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The nonlinear perturbation theory of gravitational instability is extended to include effects of both biasing and 
redshift-space distortions, which are inevitable in predicting observable quantities in galaxy surveys. Weakly 
nonlinear effects in galaxy clustering on large scales recently attracted a great interest, since the precise deter- 
mination of scales of baryon acoustic oscillations is crucial to investigate the nature of dark energy by galaxy 
surveys. We find that a local Lagrangian bias and redshift-space distortions are naturally incorporated in our 
formalism of perturbation theory with a resummation technique via the Lagrangian picture. Our formalism 
is applicable to any biasing scheme which is local in Lagrangian space, including the halo bias as a special 
case. Weakly nonlinear effects on halo clustering in redshift space are analytically given. We assume only a 
fundamental idea of the halo model: haloes form according to the extended Press-Schechter theory, and the 
spatial distributions are locally biased in Lagrangian space. There is no need for assuming the spherical col- 
lapse model to follow the dynamical evolution, which is additionally assumed in standard halo prescriptions. 
One-loop corrections to the power spectrum and correlation function of haloes in redshift space are explicitly 
derived and presented. Instead of relying on expensive numerical simulations, our approach provides an analytic 
way of investigating the weakly nonlinear effects, simultaneously including the nonlinear biasing and nonlinear 
redshift-space distortions. Nonlinearity introduces a weak scale dependence in the halo bias. The scale de- 
pendence is a smooth function in Fourier space, and the bias does not critically change the feature of baryon 
acoustic oscillations in the power spectrum. The same feature in the correlation function is less affected by 
nonlinear effects of biasing. 

PACS numbers: 98.80.-k, 95.35. +d, 95.36,+x, 98.65.-r, 



I. INTRODUCTION 

The nonlinear perturbation theory of gravitational instabil- 
ity has recently attracted renewed interest. As precision mea- 
surements of the large-scale structure of the universe become 
possible, theoretically accurate modeling is essential to inter- 
pret the observational data. It is recognized that the linear 
perturbation theory is not satisfactory for this purpose. The 
importance of nonlinear perturbation theory resides in the era 
of precision cosmology. 

A strong motivation for developing the perturbation theory 
is to model the baryon acoustic oscillations (BAOs) imprinted 
in the large-scale structure. Acoustic waves which propagate 
in the baryon-photon plasma of the early universe freeze out at 
the recombination epoch, and the sound horizon at that epoch 
is imprinted in spatial fluctuations of photons and baryons [ 1 , 
2]. The BAOs provide a standard ruler [3] to geometrically 
investigate the expansion history of the Universe fl. The dark 
energy component is efficiently constrained by galaxy surveys 
of an intermediate- to high-redshift universe (5l|g|. Using the 
BAOs as a standard ruler, large galaxy surveys are expected to 
provide a robust, promising way of constraining the nature of 
dark energy Isfj - Recent observations of BAOs in modern 
galaxy surveys [9] prove the method works well. 

Even though the BAO scale is quite large ~ lOO/^'Mpc, 
detailed structure of BAOs in galaxy clustering is affected 
by gravitationally nonlinear evolution after the recombination 



'Electronic address: taka@a.phys. nagoya-u.ac.jp 



epoch. The BAO signature in the power spectrum or in the 
correlation function is deformed by nonlinearity in a lower 
redshift universe where realistic galaxy surveys are possible 

us mm El [3. 

In galaxy redshift surveys, the clustering of dark matter is 
not directly observable. There are two major sources in the 
difference between the clustering pattern of dark matter and 
that of galaxies: galaxy biasing and redshift-space distortions. 
The spatial pattern of galaxy distribution is not necessarily the 
same as that of dark matter, and the galaxies are biased trac- 
ers of mass HH, [HI]. The redshift of a galaxy does not purely 
reflect the Hubble flow, and the Doppler shift by a peculiar 
velocity is inevitably added. Thus the spatial pattern of clus- 
tering of galaxies is distorted in redshift space lfl^[l8ll . In the 
linear regime, the power spectrum of biased objects P bj(£) 
is usually assumed to be proportional to that of mass P m (k) 
in real space: P bj(£) = b 2 P m (k), and there is a theoretical 
reason, which is known as "the local bias theorem" lfl9l r20ll . 
Assuming linear dynamics, a linear bias, and a linear velocity 
field, the power spectrum of biased objects in redshift space is 
given by IH 

Pf{k) = b 2 {\+^ i fpjk), (1) 

where b is the linear bias factor between galaxies and mass, 
P ^ Q^/b, and P m (k) is the power spectrum of mass in real 
space. 

The linear formula of Eq. ([TJ applies only when the red- 
shift is large enough, or when the wave number k is small 
enough. However, the linear formula in redshift space is not 
sufficiently accurate for typical redshift surveys 12111 . Extend- 
ing the linear formula of Eq. ([TJ to include nonlinear effects 
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is far from trivial. A straightforward application of the non- 
linear perturbation theory 1221 l23i l24l l25h with a local biasing 
scheme S has been studied so far J21 HH IHH . In the 
local biasing scheme, the number density of biased objects 
is a local function of the smoothed density field of mass in 
Eulerian space. This scheme seems to work as long as the 
tree-level perturbation theory is adopted [31]. However, when 
loop corrections are considered, there appears a conceptual 
problem in the local biasing scheme. It turns out that the one- 
loop corrections, even on large scales, strongly depend on the 
artificial smoothing scale [27]. In Ref. 12911 . a way to remove 
the dependence on the smoothing scale by renormalizing bias 
parameters is proposed. While this could be a solution to phe- 
nomenologically represent the power spectrum of biased ob- 
jects B32I1 . it is not guaranteed that the procedure actually re- 
flects the physical nature of biasing, since the biasing is more 
or less a nonlocal process. 

It is true that the biasing is difficult to be included exactly in 
an analytic framework, since the galaxy formation is a highly 
nonlinear process. The next best thing is to find a good ana- 
lytic model. In this respect, a plausible model is provided by 
a halo approach llllIllllllilllllMIMIiililll, which 
is based on the extended Press and Schechter (PS) model 
1I42I I4311 . In this approach, the galaxy biasing is described 
through two steps. In the first step, the formation and cluster- 
ing of dark matter haloes are analytically modeled by using 
the extended PS model. The spherical collapse model is com- 
bined to include the gravitational evolution of halo positions. 
In the second step, the distribution of dark matter or galaxies 
within haloes is empirically modeled with a number of as- 
sumptions [Ulimilli^liillill], which include complicated 
physics such as gas cooling, star formation, feedback effects 
from supernovae, and so forth [45]. 

While the second step involves many parameters which 
should be fitted by observations or simulations, the first step 
is theoretically less uncertain because of its purely gravita- 
tional nature. On large scales, the second step is not important 
and the power spectrum of nonlinear objects is approximately 
the same as that of haloes. Therefore, determining the power 
spectrum of haloes is an important step toward understand- 
ing the nonlinear structure formation and biasing. In the usual 
halo approach, the halo cluste ring is modeled by linear dy- 
namics and linear bias factors i33l l37l 138113911401. Recently, 
an attempt to incorporate the halo bias with nonlinear pertur- 
bation theory in a framework of local biasing scheme has been 
made [30]. There still remains conceptual arbitrariness in the 
latter formulation, mainly because the halo bias is intrinsically 
nonlocal in Eulerian space and does not fit well into the local 
Eulerian biasing scheme. From the construction, the halo bias 
is local in Lagrangian space l33ll . and is therefore nonlocal in 
Eulerian space, because positions of mass and biased objects 
are displaced by dynamical evolution. 

In this paper, a natural approach to incorporate the halo 
bias and redshift-space distortions with nonlinear perturba- 
tion theory is newly developed. To account for the locality 
of halo bias in Lagrangian space, our formulation is based 
on the Lagrangian perturbation theory (LPT) [46], instead of 
the standard Eulerian perturbation theory (EPT). In the origi- 



nal halo approach, the spherical collapse model is adopted to 
take into account the dynamical evolution in Eulerian space. 
However, the dynamical evolution of haloes is more naturally 
described in the Lagrangian picture [47]. It is not straightfor- 
ward to obtain an analytically useful formalism which com- 
bines the LPT and the halo approach. A main obstacle is that 
the observable quantities reside in Eulerian space, while the 
calculations in LPT directly give Lagrangian quantities. In 
paper I 114 811 . a new approach is developed to overcome this 
point, partially expanding the Lagrangian variables in Eule- 
rian space. The resulting expression contains an infinite series 
of perturbations in terms of the EPT, and this approach offers 
a simplified technique of resumming cosmological perturba- 
tions, such as done in renormalized perturbation theory and its 
variants ||49[ While our technique is not suitable enough 
to describe the fully nonlinear regime which other renormal- 
ization techniques aim at, we have shown in paper I that our 
technique is accurate enough in the quasilinear regime, and 
most importantly, predictions in redshift space are straightfor- 
ward. To date, our technique offers the only way of obtaining 
resummed power spectra in redshift space. In this paper, we 
show that the local Lagrangian bias, including the halo bias, 
is also straightforward to be incorporated in the approach of 
paper I, on top of the redshift-space distortions. The local La- 
grangian biasing scheme is not equivalent to the local Eulerian 
biasing scheme, and the conceptual problem about the strong 
dependence on smoothing scales which appears in the local 
Eulerian biasing scheme described above is not present in our 
Lagrangian approach. 

The rest of the paper is organized as follows. In Sec. UU 
our basic formalism is described. This formalism is applica- 
ble not only to the halo bias but also to a local Lagrangian 
bias in general. Linear and one-loop results are presented. In 
this general theory, we have parameters which are related to 
the Lagrangian biasing scheme. In Sec. [HI] those Lagrangian 
bias parameters are calculated from the fundamental concept 
of the halo model. Effects of halo bias and redshift-space dis- 
tortions on BAO scales are demonstrated both in the power 
spectrum and in the correlation function. Our conclusions are 
summarized in Sec. [IV] In the Appendix, details of one-loop 
calculations in our framework are outlined. 



II. LOCAL LAGRANGIAN BIAS 

A. Nonlinear power spectrum with a local Lagrangian bias 

In this section, we develop a method to track nonlinear evo- 
lution of the Lagrangian bias. Our method in this section is 
not restricted to the halo bias model, and is applicable to any 
bias defined by a local function of linear density field in La- 
grangian space. 

In the Lagrangian approach to track dynamical evolution of 
cosmological density fields, a set of each trajectory of a mass 
element, x(q, t), where q are initial Lagrangian coordinates, 
describes the whole property of the density field. A displace- 
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ment field ^f(q, f) is defined by 

x(q,t) = q + V(q,t), 



(2) 



and is considered as a fundamental variable of the mass den- 
sity field. Since the initial mass density field is sufficiently 
uniform, the Eulerian mass density field p m (x, t) at any given 
time t satisfies the continuity relation, 



p m (x, t)d x = p m d q, 



(3) 



where p m is the global mean density of mass. 

On the other hand, the fluid elements in which biased ob- 
jects such as haloes reside are not uniformly distributed in La- 
grangian space. Therefore, the continuity relation between the 
Eulerian density field of the biased objects, p^(x, t), and cor- 
responding Lagrangian density field, p^ b -(q), is given by 



Pobj(x,t)d 3 x = pl bi (q)d 3 q. 



objV 



(4) 



The density field p^ b j{q) represents the initial distribution of 
locations where biased objects form later. Since the formation 
of nonlinear structure is too complex to be analytically de- 
scribed from the first principle, the initial density field p^iq) 
of biased objects should be given by a good model of nonlin- 
ear structure formation, such as the halo model. 

In this paper, we assume locality of the bias in Lagrangian 
space: the Lagrangian density field p^iq) is assumed to be 
a function of a smoothed linear overdensity at the same La- 
grangian position, 



6 R (q) 
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d'q'W R (\q-q'\)6 L (q'), 



(5) 



where W R is a smoothing kernel of size R, and <5l(<?) is the (un- 
smoothed) linear overdensity. We call such biasing scheme 
a "local Lagrangian bias" in this paper. As we will show 
in the next section, the halo bias 0311 is actually a special 
case of local Lagrangian bias. Another example of the lo- 
cal Lagrangian bias is the peak bias with the approximation 
of peak-background split l3, 13 EH]. The locality of bias 
in Lagrangian space does not mean locality in Eulerian space, 
because of evolutionary effects. Therefore our biasing scheme 
does not fall into the category of the local biasing scheme in 
a usual context of EPT The local Lagrangian bias is nonlocal 
in Eulerian space. 

Thus, we introduce a Lagrangian bias function F(5) by 



Pobj(9) 



: Pobj F[6 R (q)], 



(6) 



where p bj is the comoving mean density of the biased objects, 
which is common in Lagrangian space and in Eulerian space. 
This function has the following property: 

(F(6 R )) = l. (7) 

Equation (0 is equivalent to the following equation: 

PobjW = Pobj J d 3 qF[6 R (q)] 6 3 D [x -q- ¥(g)], (8) 



where 6^ is 3-dimensional Dirac's delta function, and we sup- 
press the time-dependence for notational simplicity. Using the 
Fourier transform of this equation, we obtain an expression of 
the power spectrum of biased objects in Eulerian space, 

P obj (fc) = J d'qe~ ik M d ± d JlF(A,)F(A 2 ) 

x ^ e 'Pi5«(qi)+'i2*«(q2)]-'fc [*(qi)-*(q2)]^ _ j (9) 



where F(A) is the Fourier transform of F(6) and q - qi - qi- 
The quantity in the ensemble average (■ ■ ■ ) in the above equa- 
tion is a function of only q because of translational invariance. 
In the absence of bias, F — 1 and F(A) = 2ji5d{A), the Eq. (O 
reduces to a known expression [52]. We do not assume rota- 
tional invariance for allowing our analysis to include redshift- 
space clustering. Our convention of the power spectrum is 
given by 

(Sobj(fc)§obj(k')) = (270 3 4(fc + k')P ob j(k), (10) 



where 



^obj(fc) 



-/ 



d xe 



ikx 



Pobj 



(ID 



Similarly, the linear power spectrum Ph(k) is defined by a sim- 
ilar equation to Eq. ( fTUb and a Fourier transform of d^iq)- 

The expression of Eq. (O has a form that we can apply to 
the cumulant expansion theorem lf53ll 



( e - ix ) = exp 



s 



(12) 



where {X N ) C denotes a cumulant of a random variable X 12511 . 
The corresponding factor in Eq. (0 thus reduces to 



^ e ilM6i l (q,)+A 2 S R (q 2 )]-ik['i/(q l )-^(q 2 )]^ 

= exp 



Z 



n\+n2+m\ J rm2>\ 



n\ \n2\m\ hri2\ 



A^A^BZU^q) 



(13) 



where the multinomial theorem is used, and 



X <[<5«(gi)]"'[«q2)]" 2 [fc • *(9i)] ffll [fc • *(92)P> C ■ (14) 

The translational invariance and the parity symmetry imply 
the following identities: 

BftXft, q) = (-D m B n ^(k, -q), (15) 



= {-l) m BZLik,-q), 



(16) 



and therefore Eq. (TBI is symmetric with respect to its indices: 
(k,q) = B%"i(k,q). (17) 



D»l«2 
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When the initial density field is random Gaussian, which is as- 
sumed throughout this paper, the Equation ( fT4l > of m.\ = m.2 = 
survives only when n\ + n.2 -2: 



BT\k,q) 



&(M)> «1 = "2 = l, 



0"S> 



(o. 



(«i = 2, «2 = 0) or («i = 0, «2 = 2), 
otherwise, 

(18) 

where cr 2 = ^«(0), and is the smoothed linear correlation 
function of the linear density field, 



2^" 



j (kq)W 2 (kR)P L (k), 



(19) 



where y'o(x) = x 1 sin x is the spherical Bessel function of 
zeroth order, and 



WW) = 4n x dxj l) (kx)W R {x) 



(20) 



is the window function of the smoothing kernel. For the equa- 
tion ( TT4T > of n i = «2 = 0, we have 

^2m(fc), (wii — 2m,ni2 — 0) or (mi = 0, m2 = 2m), 
= ' B mm2 (k, q), nt\ > 1 and»2 - 1, 
0, otherwise, 

(21) 

where m is a positive integer, and 

A 2m (fc) = ([fe-*(0)] 2 '") c , (22) 
B mim2 (k, q) = (-ir {[k • *(qi)]"" [fc ■ *(q 2 )] m2 > c . (23) 

Using the above properties and quantities, and substituting 
Eq. ( fT3b into Eq. (0, we obtain an expression, 



P bj(fe) = exp 



2 Z 



(-l) m 
(2m)! 



A 2m (fc) 



af 3 #e lfeq exp 



/. j ;B mimi (k,q) 

nit }m.t 



oo 2n 2n 



F(A l )F(A 2 )e 



exp 



00 00 



■« 1 +«2 +m 1 



-Ai42&(M)+ y y — : — : ,-, ",„,„„• 



n\+tiq>\ mi+m2>l 



(2w) 3 tf3(fc). (24) 



So far the expression is formal and holds even on strongly 
nonlinear scales. Gravitational nonlinear effects on the matter 
distribution and on the bias are all included in the cumulants 
of Eq. | [1 41 . When the objects are unbiased, F(8) = 1, F(A) = 
2n6o(A), the second line of Eq. $2M simply reduces to unity, 
and the expression is equivalent to the one that was previously 
derived in Eq. (8) of paper I |48il . In the work of paper I, we 
applied the LPT to evaluate the cumulants of Eq. ( fT4b when 
«i = «2 = 0, and showed that expanding only the exponential 
factor in the integrand and keeping the exponential prefactor 
result in partial resummation of Eulerian perturbations, which 
improves the standard EPT in quasi-linear regime. This is 
justified by the fact that the exponential prefactor consists only 
of cumulants of the displacement field at a single point, while 
the remaining exponential factor consists of cumulants at two 
points separated by |q| ~ |fc| _1 . The latter cumulants are small 
enough in a large-scale limit, |fe] — > 0. When the exponential 
prefactor is expanded as well, this approach gives equivalent 
results to the standard EPT. Evaluation of the cumulants in 
redshift space is straightforward in the framework of LPT. 

In the presence of Lagrangian biasing, a similar technique 
can be adopted. In the integrand of A\ and A%, we expand 
the last exponential factor, keeping the exponential prefactor 
which involves <j 2 R . In a large-scale limit, <x 2 » £r(M) be- 
cause \q\ is large, and cr R » Bm,m 2 ( m i + m 2 S 1) because 



|fc| is small. Therefore, it is desirable to keep the exponential 
prefactor in the integrand of A's. Expanding the last exponen- 
tial factor, the integrals of A's can be performed as 



J -c 



2n 



1 r 



dde 



d"F 
~d6" 

(25) 



The right-hand side (RHS) of Eq. (125 1 corresponds to the ex- 
pectation value of derivatives of F. In the following, we use 
notations such as (F'), (F") for n — 1,2, respectively, to rep- 
resent the above integral. The local Lagrangian bias is fully 
characterized by a series of these parameters (F^), which we 
call local Lagrangian bias parameters. 

The exponential prefactor in Eq. d24T > corresponds to the 
characteristic function of the one-point distribution of dis- 
placement field. In fact, from Eqs. (fT2l) and (1221 and parity 
property, A2 m +\ = 0, we have 



exp 



2 I 



(_iy« 
(2m)! 



\(e- ik *)\ 



(26) 



Evaluation of the above characteristic function of the RHS 
requires the fully nonlinear dynamics. In this paper, the ex- 
ponent of the LHS is evaluated by adopting the perturbation 
theory. 
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Before closing this subsection, we comment on the role of 
smoothing radius R. In the local Lagrangian biasing scheme 
with peak-background split, the number density of biased ob- 
jects is spatially modulated by the linearly extrapolated back- 
ground density field in Lagrangian space. The large-scale 
clustering of biased objects should not depend on the artifi- 
cial choice of smoothing radius R to define the background 
field. In fact, in the case of halo bias, it is explicitly shown 
that (F (n) ) is independent on R. In our derivation, we have 
used an approximation cr 2 R » £r(|<zI), which means our equa- 
tions are valid on larger scales than the smoothing radius of 
the Lagrangian bias: k <K R l . In particular, the smoothing 
kernel W{kR) in final expressions of our perturbation theory is 
replaced by unity, for consistency with our approximation in 
the first place. 



B. Biased power spectrum from Lagrangian perturbation 
theory 

We evaluate the general Eq. d24t via the LPT as outlined in 
the previous section. Details of derivations by the one-loop 
perturbation theory are given in the Appendix. In this sub- 
section, we summarize basic assumptions made in the deriva- 
tions. 

As explained above, we need to evaluate the cumulants 
Bm"m 2 (k, q) of Eq. d). In LPT, the displacement field is ex- 
panded as a perturbative series, 

* = ^r(D + ^(2) + ^P) + . . . f (27) 

where is given by integrations over nth product of lin- 
ear density contrast 6l with kernels as in Eq. ( lA2t . By means 
of LPT, the cumulant of Eq. ( TT~4-b reduces to an infinite sum 
over cumulants of linear density contrast, which are straight- 
forwardly given by the linear power spectrum. On large scales 
where |fc| <K |\&|, contributions from higher-order perturba- 
tions in Eq. d27l ) are small enough, and one can truncate the 
series. In the perturbation theory, a consistent manner of trun- 
cation is given by a loop-expansion l25ll to obtain a nonlinear 
power spectrum. For Gaussian initial conditions, the loop ex- 
pansion is equivalent to the series expansion in terms of the 
linear power spectrum Ph(k). 

The treatment of nonlinear redshift-space distortions is sim- 
pler in LPT than that in EPT as shown in paper I J48ll . The 
displacement field in redshift space is given by 

z ■ 4? 

* s = * + 2, (28) 

H 

where \& is the displacement field in real space, a dot denotes 
the derivative with respect to the cosmic time t, z is a unit 
vector along the line of sight, H = a/a is, the time-dependent 
Hubble parameter, and a(t) is the scale factor. The relation 
between displacement fields in Eq. d28l i is exactly linear even 
in the nonlinear regime. In contrast, the redshift-space dis- 
tortions of Eulerian variables are given by nonlinear transfor- 
mations. This is a reason why nonlinear redshift-space distor- 
tions are easier to handle in LPT than that in EPT. 



The time dependence of each perturbative term in Eq. (1271 ) 
is approximately given by VP™ oc D", where D(t) is the lin- 
ear growth rate. This relation is exact in the Einstein-de Sit- 
ter model, and approximately holds in general cosmology 
1251, l54ll . We also apply the distant-observer approximation in 
which the line of sight z is fixed. The latter approximation is 
commonly used in analyses of redshift-space distortions and 
valid as long as the redshift surveys are deep enough so that 
clustering scales of interest is smaller than distances between 
the observer and galaxies [55]. With those approximations, 
order-by-order linear transformations of displacement fields 
become particularly simple: \I> s( "' = 7;(")^r(") ? where R <n ^ is a 
3x3 matrix whose components are given by Eq. JA7b . 

Keeping the one-loop LPT exact, the integrand of Eq. ( l24l ) 
turns out to be a strongly oscillating function of q. It seems 
extremely difficult to numerically evaluate such integral. In- 
stead, we further expand and truncate the exponential factors 
in the integrand at the one-loop level as explained in the previ- 
ous subsection. This means that our result is not exact at one- 
loop LPT, while the neglected terms are of order 0[Pl(Ic)] 3 , 
which are two- or or higher-loop contributions in terms of Eu- 
lerian perturbations. Not expanding the exponential prefactor 
improves the standard EPT as shown in paper I. 

C. The linear power spectrum in real space and in redshift 
space 

Expanding Eq. ( l24l i and keeping only linear terms in Ph(k), 
we obtain the biased power spectrum in linear perturbation 
theory. The result in real space is simply given by a linear 
term of Eq. JA691 > with a substitution / = 0, 

P obl (k) = (I + (F')f P L (k). (29) 

Since the mass power spectrum is given by P m (k) = P^(k) in 
linear theory, the linear bias factor b, which is defined by 

P obl (k) = b 2 P m {k\ (30) 

is scale-independent: 

b=l+(F'), (31) 

i.e., b does not depend on k in a large-scale limit. In the orig- 
inal halo approach, it is derived that the Eulerian linear bias 
factor is given by the Lagrangian linear bias factor plus unity 
by using a spherical collapse model j33l [34l l36ll . The result 
of Eq. OTT i. which is derived without assuming spherical col- 
lapse, is consistent to that approach. In this sense, the factor 
{F') corresponds to a Lagrangian linear bias factor. 

It is interesting to notice that the linear bias should be scale- 
independent in this limit, for any nonlinear function of F: in 
linear perturbation theory, scale dependence cannot be pro- 
duced by any form of local Lagrangian bias. This result can 
be considered as a generalization of the "local bias theorem" 
11191 12011 . which states that the linear bias factor of local Eu- 
lerian bias for sufficiently small k is scale independent. The 
constancy of the linear bias factor is now proven even for the 
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local Lagrangian bias, which is nonlocal in Eulerian space. It 
is known that the additional constant term arises from small- 
scale inaccuracies of the linear power spectrum, and the gen- 
eral asymptotic form of biased power spectrum is given by 
Pobj(k) = b 2 P m (k) + c in a large-scale limit B20I1 . 

The corresponding linear result in redshift space is given by 
a linear term of Eq. ( |A69t , 



P®£k) = (l+(F')+ffx 2 ) P L (k), 



(32) 



where fi = z ■ kfk is the direction cosine of the wavevec- 
tor k with respect to the line of sight z, f = dlnD/dlna = 
(HD) l D is the logarithmic derivative of the linear growth rate 
D(t). This result is equivalent to the Kaiser's formula 1 18] 



P®(k)=b z (l+pffp m (k), 



(33) 



where the linear bias factor b is given by Eq. OTb and fi - f/b 
is the redshift-space distortion parameter. Again, it is inter- 
esting to notice that we have derived the Kaiser's formula in 
the presence of any nonlinear local bias in Lagrangian space: 
the Kaiser's formula with a scale-independent bias is a general 
consequence of a large-scale limit even in this framework. 



D. One-loop corrections to the biased power spectrum 

The formal expression of Eq. ( t24b is evaluated by applying 
the LPT. The derivation of one-loop corrections to the power 
spectrum is detailed in the Appendix. The power spectrum in 
real space with one-loop corrections is given by putting / = 
in Eq. dA69l ). The result is 



P ob #) = exp[- (k/k NL ) 2 ] 



x |(1 + (F')fP L (k) + ^GiW + 3 -Qi(k) + ^Qi(k) 



+ {F') 



Q 5 (k) + 2Q 7 (k) 



+ {F") 



;Qs(k) + Q9(k) 



+ {F'Y [Q 9 (k) + Q u (k)] + 2(F'){F")Q u (k) 
+ \{F") 2 Qn{k) + ^ (1 + {F')f [R x {k) + R 2 (k)] 



— (i + iF'^Rmh 



where 



■tNL 



6tt 2 / 



dkP L (k) 



-1/2 



(34) 



(35) 



and the functions Q„(k), R„(k) are given by Eqs. ( IA39I )— (lA48h . 
and are second order in P^ik). When the exponential prefac- 
tor is expanded and only second order terms in P\Xk) are re- 
tained, we obtain an expression of EPT without any resum- 
mation of higher-order perturbations. In an unbiased case, 
(F') = (F") = 0, the expression reduces to the result of the 
one-loop perturbation theory of mass 1I24I l48tl . 



The power spectrum of the biased objects in redshift space 
with one-loop corrections is given in Eq. ( 1A69I >: 



J*> (fc) = exp {- [l + /(/ + 2)f?] (k/k NL f) 
X (l + (F') + tff P L (k) + Yj^TEnmik) 



(36) 



where E nm (k) is given by Eqs. ( IA7U -( lA79b . When the bias 
is not present, (F') = (F") = 0, this result reduces to the one 
derived in paper I j48ll . 

A cross power spectrum of differently biased objects is sim- 
ilarly given. When the bias functions of these objects are F\ 
and Fz, the cross power spectrum is given by substitutions of 
Eqs. ( IA18b — ( tA22l >. after expanding Eq. (f34l or d36*b in terms 
of (F') and (F"), The spherical average of Eq. (f36t can be 
obtained by using the following integral: 



= (-D" 





erf(x 1/2 ) 


[dxj 


x 1 / 2 



(37) 



where y(a, x) is the lower incomplete gamma function, and erf 
is the error function normalized by erf (+00) = 1. The correla- 
tion function is obtained by numerically Fourier transforming 
Eq. (f36b . Spherically averaged correlation function is simply 
given by 



-r 

Jo 



k 2 dk 
2^ 



Mkr)P(k), 



(38) 



where P(k) is the spherically averaged power spectrum. 

As described in paper I, the origin of the exponential prefac- 
tor exp[-(k / k^i) 2 ] is the nonlinear smearing effect by random 
motions of mass elements. In redshift space, additional smear- 
ing effect is present along the lines of sight. The latter ef- 
fect is similar to the nonlinear fingers-of-God effect 01711 . The 
form of exponential prefactor coincides with the one which 
has been phenomenologically introduced in previous work 
lfl2[ l2ltl to represent the smearing effects. 

The one-loop corrections in Eqs. d34l > and (|36| | are given by 
^-dependent functions, Q\{k), . . . , Qn{k), R\(k),R2(k), which 
are defined by Eqs. flA39MA481 l in the Appendix. All of these 
functions but Quik) vanish in a large-scale limit k — > 0. The 
function Qn{k) contributes only when (F") + 0. Thus, the 
one-loop contributions are present even in a large-scale limit 
through Q\3(k) when (F") + 0. This function turns out to be 
a convolution of the power spectrum, 

Cis(*) = J ^P^P)Pd\k-p\) = j d 3 xe- ik x [£ L «] 2 , 

(39) 

where %l(x) is the linear correlation function in real space. In 
configuration space, [£l(*)] 2 %l(x) in a large-scale limit, 
x — > 00. However, in Fourier space, Qn{k) has a finite value 
in a large-scale limit, 

Gi3(&^0) = J ^[^)] 2 = J d 3 x [^l(x)] 2 , (40) 
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while P L (k — > 0) = for cold dark matter (CDM)-like power 
spectra. Therefore, the power spectrum on very large scales 
is dominated by a constant contribution originated from non- 
linear clustering when the bias is present and (F") + 0. For 
related discussion within a framework of local Eulerian bias, 
seeRef. Q. 



E. Nonequivalence between the local Lagrangian bias and the 
local Eulerian bias 

In our local Lagrangian biasing scheme, bias parameters 
are given by a set of parameters, {(f (n) )}. In the local Eulerian 
biasing scheme, on the other hand, bias parameters are given 
by a set of parameters \b„} which are coefficients of a Taylor 
expansion, 

°° h 

*<*j = Z r"' (41) 

where 6 bj is the overdensity of objects and 6 is the evolved 
overdensity of mass at the same Eulerian position with some 
smoothing filter. One may wonder if there are some relations 
between these two sets of bias parameters. However, the two 
biasing schemes are not equivalent to each other, and the two 
sets of bias parameters are not expressible from one another in 
general. Therefore, the expression of the biased power spec- 
trum derived above is essentially different from the one with 
local Eulerian bias previously derived in literature. Below we 
clarify this situation in detail. 

An essential difference between those two schemes is that 
the Eulerian bias is applied to dynamically evolved density 
fields while the Lagrangian bias is applied to initial density 
fields. Since the dynamical evolution is generally nonlocal, 
those two local biasing schemes are not equivalent to each 
other. 

In the standard halo approach, however, the Eulerian bias 
parameters {b n } are derived l34l l4lll . although the halo bias 
falls into a category of local Lagrangian bias. One may won- 
der if the local Lagrangian bias is actually equivalent to the 
local Eulerian bias from this fact. However, the spherical col- 
lapse model is crucially assumed in such a derivation. Since 
the dynamical evolution in a spherical collapse model is lo- 
cally determined, the local Lagrangian bias and the local Eu- 
lerian bias have one-to-one correspondence in such a special 
case. It is only when the dynamical evolutions are treated 
approximately as local processes that both biasing schemes 
become equivalent to each other. 

Similarly, the linear dynamical evolution is locally deter- 
mined, and there is a relation between linear bias parameters 
of two schemes, b\ = 1 + (F'). This is the reason why the lin- 
ear power spectrum of Eq. d29l or d32l is equivalent to that 
with Eulerian linear bias. There are not such relations for 
higher-order bias parameters in general. 

Accordingly, our one-loop result cannot be obtained by just 
a reparametrization or renormalization of a set of parameters 
\b n } in one-loop EPT with local Eulerian bias. The power 
spectrum with local Eulerian bias has a strong dependence on 



an artificial smoothing length, and has a divergent result in 
a limit of small smoothing length for CDM-like power spec- 
tra rf27[ l30ll . The power spectrum with local Lagrangian bias 
derived above does not have such a strong dependence and 
one can safely take the limit. Therefore the two power spectra 
have qualitatively different properties and are never equivalent 
to each other. The strong dependence on smoothing length 
with local Eulerian bias can be removed by a renormalization 
scheme of McDonald rf29ll . Even in this case, the resulting 
power spectrum [Eq. (16) of Ref. (291] is not reachable by 
simple reparameterizations of our Eq. ( |34l , and vice versa. 
For example, McDonald's Eq. (16) has a common factor of 
b\, while the factor (1 + {F')) 2 cannot be factorized out in our 
Eq. (|34l i with any reparameterization of (F"). 

III. THE HALO BIAS IN PERTURBATION THEORY 

A. The halo approach as a local Lagrangian biasing scheme 

The nonlinear power spectrum derived above depends on 
the local Lagrangian bias only through (F ( ^) defined by 
Eq. d25l >. Up to one-loop corrections, we only need two num- 
bers, (F') and (F"). For a general local Lagrangian bias, these 
numbers could be considered as parameters which should be 
fitted by observations. Alternatively, those numbers can be 
derived once a model of bias function F(6) is specified. In 
this section, we take the latter approach, considering the halo 
bias model. 

So far the most successful biasing model in nonlinear struc- 
ture formation is provided by the halo approach JH, HH, l35l 
[^[3^E|j39l,[40il4lll, which is based on the extended PS the- 
ory 11421 14311 . In the PS formalism, the mass of halo is related 
to the Lagrangian radius R of spherical cell by M — AnpR i /3, 
or/? = [M/(1.162 x \0 n h- l M Q Q. m )] ll3 h- l Mvc, where M Q = 
1.989 x 10 30 kg is the solar mass, and Q m is the density pa- 
rameter at the present time. The variance of mass overdensity 
in the cell, as a function of mass scale, which is linearly ex- 
trapolated to the present time z — 0, is given by 

/k 2 dk 
— W 2 (kR)P (k), (42) 

where W(x) = 3(sinx-xcos x)/x 3 is the top-hat window func- 
tion, Po(k) = P^(k)/D 2 is the linear power spectrum extrapo- 
lated to the present time. 

The critical overdensity, which is required for spherical col- 
lapse at redshift z, and is linearly extrapolated to the present 
time, is given by 

£ c (0) 

SM = W> m 

where D(z) is the linear growth rate as a function of red- 
shift z, normalized as D(0) = 1 . The critical overdensity at 
the present time, <5 C (0), is given by 6 C (0) = 3(3tt/2) 2/3 /5 = 
1.68647 in the Einstein-de Sitter model. In general cosmol- 
ogy, 6 C depends weakly on cosmological parameters and red- 
shift 15611 . It is still a good approximation to use the constant 
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value of the Einstein-de Sitter universe in general cosmologi- 
cal models. 

According to the PS theory, the comoving number density 
of haloes with mass between M and M + dM, identified at 
redshift z, is given by 



2p dv 
n(M,z)dM= —g(v)—, 



(44) 



where v = S c (z)/(r(M), g(v) = (27rr 1/2 vexp(-v 2 /2), and p 
is the comoving mean density of mass. In literature, a multi- 
plicity function f(v) defined by g(v) - vf(v)/2 and v = v 2 is 
frequently introduced. The PS mass function is improved by 
Sheth and Tormen (ST) ll36ll to give a better fit to that of haloes 
in numerical simulations of CDM-type cosmologies. The ST 
mass function is also given by Eq. (l44l with a modified func- 
tion 



g(y) 



A(p) 

V2^ 



1 



(qv 2 )P 



yfqve 



-1v 2 /2 



(45) 



where A(p) = [1 + 7T 1/2 2-*T(l/2 - p)]~\ and p = 0.3 and 
q = 0.707 are numerically fitted parameters. The PS mass 
function is also given by Eq. (f43T > with p = and q — 1 . The 
Eq. d45l ) has the following normalization, 



dv 

g(y)— 

o v 



(46) 



which is equivalent to \n(M, z)MdM = p. 

The original PS theory is extended to give the number den- 
sity of haloes of mass M\ , identified at redshift Z\-, in a region 
of Lagrangian radius Ro in which the linear overdensity ex- 
trapolated to the present time is So'- 



n(M l ,z l \6 ,Ro)dM l = ^.^10)—, 
Mi vio 



where 



vio = 



61 - 60 



.1/2' 



K-0 

61 = 6 c (zi), 0-0 = o-(Mq), a- 1 = o-(Mi) 



(47) 

(48) 
(49) 



and Mo = 4-7rpR^/3. The haloes of mass Mj are collapsed at 
Zi, while Mo is assumed uncollapsed at z = 0, and thus we 
always have S\ > So- 

The conditional number density of Eq. (J47j is interpreted 
as a large-scale spatial modulation of halo densities in La- 
grangian space. The number density of haloes is a function 
of the linearly extrapolated overdensity with smoothing ra- 
dius Rq. Therefore, Eq. d47l > corresponds to the Lagrangian 
number density of biased objects in Eq. ((6) for haloes of mass 
between Mi and Mi + dM\ . The bias function for the halo 
bias is given by 



F Ml (So) 



n{M u z\\So,Ro)dM x 
n(M u zi)dM\ 



(50) 



where the mass scale Mi is explicitly denoted. Defining v\ = 
<5i/o"i, we have 



dvio/no 
dvi/vj 



en 



2 2 ' 



(51) 



where v\, v\o, <J\, cr are considered as functions of halo mass 
Mi. Combining Eqs. d44l . WIS . ( l50b . and (TJTJ, the halo bias 
function reduces to 



F Ml (So) 



g(vw) 



(52) 



Since the RHS of Eq. 
of vio, we have 



07 -o-q y (vi "> 

depends on 60 only through Si - So 



(53) 



where the average (■ ■ ■ ) is taken over distribution of So- As- 
suming the initial density field is Gaussian, the distribution of 
So is also Gaussian with variance <x 2 . In the case of PS mass 
function, p — and q — 1, a Gaussian integral to give the 
average {g(v\o)) can be exactly performed: 



<0(vio)> 



1 r°° 

——= d6oe- s ° /2cr °g(v m ) 



2 2 



g(v\). 



(54) 



In a general case of ST mass function, p + and q + 1, 
Eq. (l54l i is shown in a limit cr <k en, which is a reasonable 
approximation when Mo is much larger than M\ . 
From Eqs. (|52]|-(|54]|. we have 



(-D" / 9 
g(vi) \dSi 



(-\y Vi n d n g{vi) 

9(n) = cm/ \ n ' ( 55 ) 
(5l gr(vi) dvi" 



With this expression, the following consistency relation is 
straightforwardly shown: 



1 r°° 

-J dM x n(M u zi)M x (F^) = 




(56) 



For our purpose of one-loop corrections, we need only first 
two derivatives. Dropping the subscript M and '1', they are 
given by 



(F') = 



(F") = 



1 



Sc(z) 

1 

Sj(z) 



qv 



1 + 



2p 



3qv' 



(qyiy 
2p (2qv 2 +2p - l) 



1 + (qv 2 ) 



2\l> 



(57) 
(58) 



where v = S c (z)/cr(M). The Eulerian bias factor of Eq. (f3Tb 
with Eq. d57l i agrees with that of the original halo approach 
II3HI3HI36I1. In Fig.ffl {F') and {F") are plotted against mass 
of haloes. We adopt cosmological parameters Q m = 0.28, 
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M [M Q /h] 



FIG. 1: Local Lagrangian bias parameters (f), {F") as functions 
of halo mass. Different curves correspond to different redshifts (z = 
0, 0.5, 1, 3 from bottom to top in each panel). 



O a = 0.72, D. h = 0.046, h = 0.7, n s = 0.96 and cr 8 = 0.82 
and the linear power spectrum Po(k) is calculated from the 
output of the camb code B57I1 . This set of parameters is always 
assumed in the following figures throughout this paper. 

When a finite range of mass [M\,M2\ of haloes is consid- 
ered, the denominator and the numerator of Eq. ( T50b are al- 
tered into integrations in that range. The above derivation is 
similarly applied in this case. As a result, we have 

f M2 „ d"gd In cr dM 

i m\ (-1)" Jm, dv n dM M 

\ / 6«(z) C M > , diner dM ' ' 

I o(v) 

J M| yW dM M 

In a limit M2 — > Mi, the Eq. (l55l l is recovered as expected. 



B. Scale-dependence of the bias in quasi-linear regime 

Applying Eqs. ( l57l i and ( f58l ) to our results of Eq. d34b or 
Eq. (|36l l, a power spectrum of haloes with one-loop correc- 
tions, either in real space or in redshift space, can be evaluated. 
To demonstrate the effects of biasing and redshift-space dis- 
tortions, we show examples of power spectra and correlation 
functions in this subsection, and briefly discuss the impact on 
BAO features. Since the main purpose of the present paper is 
to give a new formalism in perturbation theory, detailed inves- 
tigation of BAOs with our approach will be given elsewhere. 

In Fig. |2] the normalized power spectra in real space are 
plotted, and those in redshift space are shown in Fig. [3] We 
adopt the same set of cosmological parameters as in Fig. [TJ 
Angular averages are taken for power spectra in redshift 



space. In paper I, it is shown that our one-loop formula is 
valid within a few percent for k < &nl/2 B48I1 . compared to 
numerical simulations. Expecting this criterion is also effec- 
tive in our generalization including bias, the corresponding 
scales of validity, &nl/2, are shown in vertical dotted lines 
in the figures. To highlight nonlinear effects, overall ampli- 
tudes predicted from linear growth rate, linear bias and lin- 
ear redshift-space distortions are scaled out. In real space, 
the amplitude of power spectrum is proportional to a scale- 
independent factor D 2 (z)b 2 (z), where b(z) is the linear bias 
factor defined by Eq. ([3TT >. In redshift space, Kaiser's en- 
hancement factor R(z) = 1 + 2/3/3 + y3 2 /5 iQji is an additional 
source of the linear amplitude. 

In upper panels in the figures, each power spectrum is nor- 
malized by a smoothed, no-wiggle linear power spectrum 
P nvl (k) of Ref. B58I1 to highlight baryonic features. The lin- 
ear amplification factors described above are all scaled out. In 
lower panels, power spectra of haloes divided by those of dark 
matter are plotted, where amplifications by linear biases and 
linear redshift-space distortions are scaled out. Thus, curves 
in lower panels show the scale dependence of bias. In the 
usual halo approach, the scale dependence of the bias arises 
only from galaxy /dark matter clustering within haloes, which 
is not considered in this paper. The scale dependence shown 
in our results purely originates in clustering of haloes them- 
selves. Linear theory predicts constancy of halo bias on large 
scales. Nonlinear effects of dynamics, biasing, and redshift- 
space distortions are responsible for the scale dependence. 

Comparing the clustering in real space and in redshift 
space, the power spectra on small scales are suppressed by 
nonlinear redshift-space distortions. This suppression is due 
to the large-scale random motion of objects, which is simi- 
lar to a phenomenon known as a fingers-of-God effect 11711 on 
small scales. Generally, the scale dependence of bias is strong 
for very massive and very light haloes. Haloes of intermediate 
mass do not show significant deviations from constant bias. 
The scale dependence of bias does not show significant oscil- 
lations, and are mostly smooth functions of scales. This is de- 
sirable for cosmological applications to use the BAO scale as 
a standard ruler. Various nonlinear effects mostly modify the 
broadband shape of power spectrum, and therefore resulting 
shifts of acoustic scales are correctable as numerically demon- 
strated in Ref. JH. 

Although the power spectrum and correlation function are 
related by Fourier transforms and have mathematically equiv- 
alent information, cosmological information that can be ex- 
tracted from them with real data is not exactly equivalent to 
each other, because error properties are different. While the 
BAO scales are imprinted in multiple wiggles in the power 
spectrum, there appears one single peak in the correlation 
function H]. The physical BAO scale is just a single scale, 
and many wiggles in the power spectrum are overtones of 
the fundamental scale of BAOs. In Figs. |4]and|5] correlation 
functions are plotted. They are calculated by Fourier trans- 
forming the power spectra of Figs. [2] and [3] In paper I, we 
have shown that our method of one-loop corrections to the 
clustering of dark matter, both in real space and in redshift 
space, agrees very well with numerical simulations on scales 
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FIG. 2: Dependencies on halo mass and redshift of nonlinear power spectrum in real space. In the top panels, each power spectrum is divided 
by a smoothed, no-wiggle linear power spectrum P nw (fe) (Hj], and by a squared linear bias factor b 2 . Values of redshifts and halo masses are 
shown in each panel. Solid lines: nonlinear power spectra of haloes with different masses with increasing order from thinner to thicker lines; 
dotted lines: linear theory; dashed lines: nonlinear power spectra of dark matter. In the bottom panels, halo power spectra are divided by 
corresponding mass power spectra and by squared linear bias factor, presenting the scale dependence of halo bias. Vertical short-dashed lines 
correspond to the scale k NL /2 to indicate the validity range k < £ NL /2, where our result is expected to be accurate within a few percent. 




FIG. 3: Same as Fig. [2] but in redshift space. Spherically averaged power spectra are plotted. Linear redshift-space enhancement factor 
R = 1 + 2/3/3 + /3 2 /5 is also scaled out. 
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FIG. 4: Dependencies on halo mass and redshift of nonlinear correlation function in real space. Correlation functions with a fixed redshift and 
with different halo masses are presented in each column. Mass of the halo varies in increasing order from thinner to thicker solid lines. Dotted 
lines correspond to the prediction of linear theory and dashed lines correspond to nonlinear correlation functions of dark matter. In the top 
rows, the bare values of correlation function are plotted. In the middle rows, the correlation functions are normalized by linear bias factors and 
linear growth factors. In the bottom rows, residual values in the normalized correlation function of haloes (plotted in middle rows), relative to 
that of dark matter, are plotted. 




M/(10"M o /h) 

= 1, 10, 10 2 , 10 s 
I | I I I | I I I | I I I 



M/(10"M s /h) 

= 0.5, 5, 50, 500 

| I I I | I I I | I I I 




M/(10"M s /h) 
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r [Mpc/h] r [Mpo/h] r [Mpc/h] 



FIG. 5: Same as Fig. [4] but in redshift space. Spherically averaged correlation functions are plotted. Linear redshift-space enhancement factor 
R is also scaled out. 
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r > 10 h Mpc, where BAO signatures appear. Upper panels 
show bare correlation functions, and middle panels show nor- 
malized correlation functions. The overall amplitudes of the 
linear growth rate, linear bias and linear redshift-space distor- 
tions are scaled out in the normalized plots. The normalized 
correlation functions do not significantly depend on mass of 
haloes, and plotted lines are quite overlapped. In lower panels, 
residual values relative to normalized correlation functions of 
dark matter are plotted. 

Nonlinear effects degrade the signature of the BAO peak 
because of random displacements of haloes, and the degrada- 
tion is larger in redshift space because of additional displace- 
ments along the lines of sight. Halo bias does not signifi- 
cantly change the shape of the BAO peak. This property is 
consistent with a recent analysis of numerical simulations of 
halo clustering [14]. Effects of nonlinear dynamics and non- 
linear redshift-space distortions dominate those of nonlinear 
halo bias. Massive haloes slightly enhance the BAO signature 
by 5% or so, while less massive haloes only change by 1-2%. 
The effects are slightly larger in redshift space. 



IV. SUMMARY 

In this paper, we show that the nonlinear perturbation the- 
ory via the Lagrangian picture is naturally incorporated with 
the local Lagrangian bias on top of the redshift-space distor- 
tions. Explicit results of one-loop power spectrum with a halo 
bias both in real space and in redshift space are given. Our 
approach does not suffer from a conceptual problem which 
appears in the EPT with local Eulerian bias. The halo bias, as 
a local Lagrangian bias, is properly treated in the Lagrangian 
picture of perturbations. Simultaneous inclusion of redshift- 
space distortions is also natural in the Lagrangian picture. Our 
approach involves a partial resummation of Eulerian pertur- 
bations, and provides a better description in the quasilinear 
regime than the standard EPT. 

In our general framework of local Lagrangian bias, the bias 
is characterized by local Lagrangian bias parameters 
We do not need a Taylor expansion of the bias function as fre- 
quently adopted in local Eulerian biasing scheme. Only the 
first two of the local Lagrangian bias parameters are needed 
in our one-loop calculations. The results in Sec. [ll]are appli- 
cable for all biasing schemes as long as the bias is local in 
Lagrangian space. The local Lagrangian bias parameters are 
considered as free parameters if the bias model is not spec- 
ified. Because of the nonlocal nature of gravitational evolu- 
tion, the local Lagrangian biasing scheme is not equivalent to 
the local Eulerian biasing scheme. 

There is a successful model of Lagrangian bias, i.e., the 
halo bias. The general halo approach consists of several pieces 
of assumptions. We adopt only the most fundamental con- 
cept of the halo approach, the halo bias in Lagrangian space. 
The clustering of galaxies or dark matter is dominated by halo 
clustering on large scales, where the perturbation theory is 
safely applied. As a result, the local Lagrangian bias parame- 
ters are calculated for halo bias without ambiguity. The result- 
ing power spectrum does not have any free parameters once 



the mass or mass range of haloes is specified. 

In usual linear analysis, the halo bias is independent on 
scales. However, nonlinear effects introduce the scale depen- 
dence into the halo bias. Such scale dependence could af- 
fect the determination of the BAO scales unless the effect is 
properly quantified. Such effects are quantitatively calculated 
without numerical simulations. We find that the scale depen- 
dence of bias is a smoothly varying function. Therefore, the 
BAO scales are shifted by change of the broadband shape in 
the power spectrum. The correlation function is also affected 
by scale dependence of nonlinear halo bias. Compared to the 
power spectrum, the shape of the BAO peak in the correla- 
tion function seems to be less affected. Those observations 
are consistent with the recent results of numerical simulations 

MM. 

Our formalism is compatible with any biasing scheme 
which is local in Lagrangian space. For example, the peak bias 
is approximately considered as a local in Lagrangian space in 
a limit of the peak-background split lfT6l [3411 . when the con- 
straints imposed by the spatial derivatives to define peaks can 
be neglected. However, the last approximation may not be 
appropriate for accurately predicting the BAO signature M60I1 . 
The exact biasing mechanism in the real world is definitely not 
local both in Eulerian space and in Lagrangian space. The suc- 
cess of the halo approach indicates the local Lagrangian bias 
is a good approximation at least on large scales. However, ex- 
tending the model of local Lagrangian bias to a nonlocal one 
is an option to make the theory more accurate and general. It 
is straightforward to extend our formalism in Sec. |II]to include 
a nonlocal bias. 

Although our formalism contains a partial resummation of 
higher-order Eulerian perturbations, Lagrangian perturbations 
are truncated at the one-loop level. Recent developments 
of the renormalized perturbation theory and its variants ll49ll 
show that it is possible to reorganize and resum higher-order 
perturbations using the concept of propagators. It would be 
interesting if one could use the concept of propagators in La- 
grangian space 116 ill to further renormalize the present formal- 
ism and to describe the deeply nonlinear regime, k > £nl/2, 
which is not accessible by the present formalism. 

We consider only the halo bias in this paper. The result- 
ing power spectrum corresponds to that of halo centers, or the 
two-halo term of the halo model on large scales. In the halo 
approach, the nonlinear power spectrum is given by a super- 
position of the one-halo term and the two-halo term, which 
are given by convolutions with a model of density profile of 
galaxies or dark matter. A model of density profile is domi- 
nantly relevant to clustering on small scales. In a context of 
the halo approach, our formalism improves the description of 
the two-halo term, in which only the linear dynamics is usu- 
ally included because of simplicity. It is possible to include a 
model of density profile of halo model in redshift space B62I1 . 

As pointed out in paper I, our approach does not have much 
power on small scales in the power spectrum, k > k^/2. The 
main source of the powerless is the exponential damping fac- 
tor, which originates from random motion of the displacement 
field. A fully nonlinear description of this factor may dra- 
matically improve the applicability of the present formalism 
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for dark matter clustering. However, in the presence of bias, 
the nonlinear regime is dominated by the scale dependence of 
bias, which may be more appropriately described by halo ap- 
proach with the one-halo term, including a model of density 
profile or halo occupation dynamics. 
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APPENDIX A: ONE-LOOP CORRECTIONS TO THE 
BIASED POWER SPECTRUM VIA THE LAGRANGIAN 
PERTURBATION THEORY 

In this Appendix, we outline a derivation of one-loop cor- 
rections to the power spectrum with the local Lagrangian bias. 
Our goal is to evaluate the cumulants of Eq. ( TPfl i by the pertur- 
bation theory, and obtain a perturbative expansion of Eq. ( 1241 ). 
Most of the necessary techniques have been developed in pa- 
per I H. 

For the one-loop corrections, the cumulants of Eq. ( fT4b up 
to second order in Piik) should be evaluated. We have applied 
the LPT [46, 48] to evaluate the similar cumulants in paper I, 
and we can use the same method. In the LPT, the displacement 
field \J/ is expanded by a perturbative series: 

* = ^r(D + *(2) + ^r(3) + . . . < ( A1 ) 

The first-order term corresponds to the classic Zel'dovich 
approximation [63]. The spatial derivatives of each term, 
di^j have the order of (9(<Sl)", where <Sl is the linear den- 
sity field. In a Fourier representation, 

nlj (2nf (2nf K (p 

x L ( "\pu . . .,p„)6 L (Pi) ■ • -<5l(P„), (A2) 

where <5l(p) is the Fourier transform of the linear density field, 
and perturbative kernels i (n) are given by the LPT. Since the 
dependence of these kernels on time and on cosmological pa- 
rameters is weak, it is a good approximation to use the ker- 
nels of the Einstein-de Sitter model even in general cosmol- 
ogy Oil]. In real space, expressions of the kernels L (n) up to 



third order are given by j48ll64ll 

= |, (A3) 

^w^—fi-^y, (A4) 

i (3a) (Pl,P2,P3)=~(l-^, 2 )(l-^ 2 ,3) 

- - 3j"i,2 + 2 / < i,a u 2,3A'3,i) + fexT(p 1 ,p 2 ,p3), 

(A5) 

where fe = pi + ■ ■ ■ + p„ for each 2/"\ fitj = p t ■ Pjlipipj), 
P-ijM - (Pi + Pj) • Pk/(\Pi + Pj\Pk), and a vector T represents a 
transverse part whose expression is not needed in the follow- 
ing application. It is useful to symmetrize the kernel L^ 3a) in 
terms of their arguments: 

L (y> (p u p 2 ,Pi) = i [i (3a, (pi,p 2 ,P3) + perm.] . (A6) 

As shown in paper I, the perturbative kernels in redshift 
space are simply given by linear transformations by redshift- 
space distortion tensors R^"\ whose components are 

B$ = Sij + nfziZj, (A7) 

where / = dlnD/dlna = (HD) l D is the logarithmic deriva- 
tive of linear growth rate D(t) by the scale factor a( t), and z, is 
a unit vector along the line of sight. In a matrix notation, the 
kernels in redshift space is given by 

L m _ flOO^M (A8) 

It is useful to define the following mixed polyspectra of lin- 
ear density field and displacement field: 

= (2n) 3 6 3 D (fei + • • • + fe; + p\ + ■ ■ ■ + p,„) 

x (-0'"C/, ..,,„ (fci , . . . , kf, p u . . . , p m ) , (A9) 

where (5l and "P; are the Fourier transforms of the linear den- 
sity field and the displacement field, respectively. When / = 0, 
the above polyspectra are equivalent to the ones defined in 
Eq. (1 1) of paper I, but we adopt an opposite sign in this pa- 
per. For I + m — 2 in the Eq. ( |A9t , we also use notations such 
as 

C(fc) = C(fe, -fe), CKfc) = C,(fc; -fe), Cy(fc) = C i7 (fe, -fe). 

(A10) 

When m = 0, the above polyspectra of Eq. ( IA9b is nonzero 
only when I = 2 for a Gaussian initial condition, which is 
assumed throughout this paper. The Eq. <fT4T > has the order 
Q\P] .(k)] ll + ' 2 +m ' + '" 2 ~ 1 , because of the property of cumulants 
112511 . Therefore, we only need to consider n\ + n 2 + m \ + 
ni2 < 3 up to one-loop corrections. Expanding the exponential 
factors in Eq. d24"i i. but the first prefactor, and truncating third- 
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or higher-order terms in Ph(k), we obtain 



P obi (k) = exp 



obj 
X 

where 



1 J (2ti) 3 ' 
[floo(fc) + <^'> flio(fc) + (F") aoi(fc) + (F') 2 a 20 (k) 

+{F')(F")an(k) + (F") 2 ao 2 (k)], (All) 



floo(fc) = - kikjCij(k) - kikjk k J ~n^f Ci ^ k ' ~ p > P~ k "> 

+ \k i k j k k k l f ^PCij{p)C u {k - p), (A12) 
2 J (Iny 



flio(fc) = 2kjCi(k) 
+ 



J [c,y(fc; -p, p-k)- 2Cij(-p; p - k, k)j 
-2kikjk k J ^Q(p)C jk (k-p), (A13) 

floi(fe) = -ki f -^-^C,(-p,p- fc;fc) 

+ k i k j f jLE-C,(p)Cj(k-p), (A14) 

a 20 (k) = C(k) + 2ki f ^-£-Ci(k,-p;p-k) 

+ ktkj J ^ [Ci(p)Cj(k - p) - C{p)C u {k - p)] , 



J (2w> 



:C(p)C,(fe-p), 



(A15) 



(A16) 



aoi(k) = l f AlLc(p)C(k-p). (All) 

We have neglected effects of the smoothing kernel W(kR) in 
the above equation, for the consistency of our treatment as 
discussed in the end of Sec. IIII Al A cross power spectrum 
of differently biased objects is similarly given. When the bias 
functions of these objects are F\ and F 2 , the cross power spec- 
trum is given by substitutions 



{F')^ l -{{F[) + {F' 2 )), 

(F")^\{(F'{) + (F'{)), 

(F'f -» (F\)(F' 2 ), 
(F"f -» {F'(){F'{), 

(f'xf") -» I «f;><f^'> + <f;'><f^» , 



(A18) 

(A19) 

(A20) 
(A21) 

(A22) 



in Eq. ( IA1U . 

Next we define mixed polyspectra of each order in pertur- 
bations: 

(&.(fcl) ■ • • ^(fc^'fo) • ■ ■ *£" W) c 

= (27T) 3 5 rj (fej + • • • + fe; + Pi + • • • + p m ) 

x (-fTCj::^ (fei, . . ., fe/;pi, . . . ,p m ) , (A23) 

where *Pr are the Fourier transforms of the displacement field 
of order n in Eq. (lAlb . For I + m — 2 in the Eq. dA23l ). we also 
use notations similar to those in Eq. ( 1A10I ). such as C (n) (fc), 
Cf°(fc), C^'" 2) (fc). The Eq. (TA231 is nonzero only when I + 
ri\+- ■ -+n m is an even number. The original mixed polyspectra 
of Eq. dA91 > are given by sums of the polyspectra of each order. 
In particular, 



dip) = Cf ) (p) + C, w (p) + 



<3), 



(A24) 



djip) = c^»(p) + cf;\p) + cf/\ P ) + cf}\p) + ■■■, 

(A25) 



Ci(pi,p 2 ;p3) = Cf'(pi,p2;p3) + 



(A26) 



Cij(PuP2,P3) = C, ( ' 2) (Pi;P2,P3) + Cf'(Pi;P2,P3) + 



(21) c 



C yi (pi,p 2 ,p 3 ) = C ( i 1 j 1 k 2) (p U P2,P3) 



(A27) 



+ cg^tpi.^j*) + Cg 1, (p 1 ,p 2 ,p 3 ) + ■ ■ • , (A28) 

up to second order in Ph(k). Using the LPT kernels of 
Eq. dA2b . the mixed polyspectra of each order are given by 



C(p) = P L (p), 

c, (1) (p) = Lf\ P )pap\ 



(A29) 
(A30) 

C?j l) (p) = -L^{p)L ( J\p)Pdp\ (A31) 



Cf ) (p) = ^P h (p) J ^Lf\p,-p\p')P^p\ (A32) 

c<f (p) = -i J ^4 2) (p'.p-p')^ 2) (p',p-p') 

xPL^'^Ldp-p'l), (A33) 

C<] 3) (p) = Cj»(p) 
= -iLP(p)P L (p) J ^Lf(p,-p',p')P L (p'), (A34) 



C, (2) (pi,p 2 ;P3) = Lf\puP2)PL(p l )PL(p2), 

C-] 2) (pi;p2,P3) = C (2 i 1 \pi;p3,p2) 



(A35) 



= -L) i >(p 2 )l};>(p 1 ,p2)P L (p l )P L (p2), (A36) 

C^ 12) (P1,P2,P3) = Cg 1) (p 3 ,Pl,P2) = C^ 1) (p 2 ,P3,Pl) 

= L? ) (p l )Lf ) (p 2 )Lf ) (p u p 2 )P L (p l )P h (p 2 ). (A37) 

As in paper I, diagrammatic representations are helpful to 
understand the structure of perturbative terms. With Feynman 
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P 



x Pi 



Li(pi,-.-,p„ 




Pi 



K---K,{F (n) ) 



FIG. 6: Feynman rules for diagrammatic representations. 



C(p): 



4 3) (P) : 



FIG. 8: Feynman rules for external lines. The momentum conserva- 
tion, k = Pi + ■ ■ ■ p m +p\ + • • • p' n , is assumed. 



P i 



Cf>(p) : 



i P .1 



C« 3, (P) 




Cfdn.pa,-!*) : 
- • Pi 

. ' p 

cf '(Pi;P2,p 3 ) : 
* ^ Pi 

— ■ — r- " Pi 

i 

-,(112) 



C.jfc (Pl,P2,Ps) 



Pi 



P2 



of the nonlinear power spectrum in quasi-linear regime than 
the standard EPT. In the diagrammatic interpretation, the use 
of exponential prefactor corresponds to the partial renormal- 
ization of the vacuum graphs, as depicted in Fig.fTUl 

In real space, it is straightforward to evaluate Eq. (IA1 II ) 
using Eqs. (IA12I> — dA37t and perturbative kernels in real space, 
L^"\ The above expressions are not confined in real space, 
and they are also applicable for redshift-space clustering when 
we use perturbative kernels X s(n) in Eqs. ( I A29b — ( I A37t . As a 
result, the mixed polyspectra in redshift space are obtained 
from those in real space, subjecting to linear transformations: 

Q(n v -n„) _^ QS(n v ~n m ) _ _ _ , J^Tht) "%) (A38) 

ivhn 'V"4» hjl imjm jvjm ' 

Therefore, it is desirable to first perform the integra- 
tions which appear in Eqs. ( IA12I >— dA17| > with decomposed 
polyspectra of Eqs. ( IA24b — (IA27t in real space. The calcula- 
tions are similar to those presented in Appendix A of paper I. 
Several integrations we need are already given there, and oth- 
ers are not. To present the results, we first define the following 
integrals: 



FIG. 7: Diagrammatic representations of polyspectra. 



rules of Fig.|6]and appropriate statistical factors, Eqs. dA291 >— 
( IA37t are diagrammatically represented in Fig. [7] Substitut- 
ing Eqs. (TaT21 -( TaT71 > and GM4\ - (A37\ into Eq. (IaTTI . we 
obtain a lengthy expression of P \,j(k). For a diagrammatic 
representation of the result, we introduce additional Feynman 
rules for external lines in Fig. [8] All the contributions to the 
power spectrum P \,j(k), but the exponential prefactor, are di- 
agrammatically given in Fig. [9] In the exponent of the expo- 
nential prefactor in Eq. ( IA1 lb . we consider only the first order 
in P\Sk), since the remaining factor is already first order. In 
this approximation, the exponential factor corresponds to the 
bubble diagrams in Fig. [10] If we expand this exponential fac- 
tor as well, we obtain a standard loop expansion of the power 
spectrum. In the absence of bias, the result agrees with that of 
standard EPT. In paper I, we pointed out that keeping the ex- 
ponential prefactor unexpanded provides a better description 



k 3 r°° 

Qn(k) = ^J drP h (kr) 



x I dxP L [k(l + r 2 -2rx) 1/2 ] Q„(r,x), (A39) 



k 3 r°° 

**<*> = ^» J drPdkr)R n (r), 



where 



(1 + r 2 - 2rx) 



x 2 (l - rx) 2 
Q3 = (l + ri-2rxr QA 



(1 + r 2 - 2rx) 2 



1 -x 2 
(1 + r 2 - 2rx) 2 



1 + r 2 — 2rx 



1 + r 2 — 2rx 



(A40) 



r 2 (\-x 2 ) 2 {l-x 2 )rx{\-rx) 
- 7, . 7 ^ — v>< = — - — ^ — — » (A41) 



, (A42) 



Qs = P^, fe ^ 1 : 3 ^ 1 :^ , (A43) 






FIG. 10: Diagrammatic representation of the exponential prefactor 
in one-loop approximation. 



27 = 



x 2 (l -rx) 
I + r 2 — 2rx 



r 2 (l -x z ) 
1 + r 2 — 2rx 



- rx(l - rx) - 2 
Q ^ l + r 2 -2rx 2k. = 1-x, 



Gll=*% Qn = rx, Qn^r 2 , 



(A44) 



(A45) 



(A46) 



and 



f r 2 (l-j 
R\ = dx- — 

J-i l+r 2 - 



x 2 f_ 
2rx 



1 + r 2 

Ha? 



(r 2 -D 4 
16r 3 



-(3 - 14r 2 + 3r 4 ) + 



In 



1 +r 



l-r 



(A47) 



Ri 



x )rx(l - rx) 
+ r 2 — 2rx 



-(3 - 2r + 3r ) + — — ; In 



24r 2 



16r 3 



1 +r 



l-r 
(A48) 



(A52) 



(A53) 



C ( ] 3) (fc) = C ( ; 31) (fc) = -l^-R.ik), (A54) 



" P Cf ) (fe,-p;p-fe) = ^[A' l (/v) +A' 2 (/v)|. (A55) 



(2tt) 



3 fc; 



Cf >(-p, p-k;k) = -~Q 8 (k), (A56) 



3 5u 3 

+ — -^[*i(*) + 2R 2 (*)], (A57) 



d 3 p not 3 fc/fc; 

-JLcf\-p;p-k,k) = ---^-Q 5 (k), (ASS) 



(2n) 



After lengthy algebra (see also paper I), we obtain 

C(fc) = PL(Jt), (A49) 



d 3 P r (2i) 

(2ff) 3 y 



C-; ^p-p-k,k) = -^[/?iW+/? 2 W], (A59) 



J (2tt) 3 * ^ ^ J (2;r 



T )3 yt 
3 kjdjk 
14 A: 4 



<^ 3 P ,-,(211) 



d 3 p_ r (i2i) 



C\f\k,-p,p-k) 



3 fc-fcfat 

Ri(k) - — -^mk) + 2R 2 (k)], (A60) 



C ;"(fe,-p,p-fe) 



_ 3 fc;<fy 

~ 14~ 



3 Art 

GiW-Tr-^teiW + ZffaWL (A61) 



cPp 



14 A: 6 



C(p)C(fc-p) = Qi 3 (*), 



^Ac(p)c ( a) (fe-p) = |e 12 (fe), 



(2tt) 3 



(A62) 



(A63) 



(2^) 3 ^ ,J ^ 

~Qm(k) + [fiio(t) - 2Qu(*)], (A64) 



2 fc 2 



2 A 4 



| ^c;"(p) C y-*-P) 



I Si 



1 



= -r 7tC2b(*) + = + 2Q 9 (fe)], (A65) 



2/t 2 



2 fc 4 



(2jt) 3 



1 5(nkk\ 1 kjkjki 

) -Q6(k) + -—±-[Q 6 (k)-2Q 7 (k)], (A66) 



2 fc 4 



2 A: 6 



i!£-C^(p)d">(fe-p) = ^^ 
(2?r) 3 ^ ^ 8 A 4 

1 S(jjkkki) 



Qi(k) 



4 £ 6 
1 kjkjk k ki 

8 /fc 8 



[3fii(*) + 12f22(t)-2(24(t)] 

[3fii(t) + 24g 2 (£) + 8g 3 (fc) - 4<2 4 (&)L 

(A67) 



where the spatial indices are symmetrized over in Eqs. (IA66b 
and ( IA67b . In deriving the above equations, the transverse part 
of Eq. ( IA5I ) does not contribute at all, because the rotational 
covariance implies 



d 3 p 



(2n) 



3 g(k, p)T(k, -p, p) oc k, (A68) 



where g is a scalar function. 

Applying the transformation of Eq. (|A38t , the correspond- 
ing integrals in redshift space are straightforwardly obtained. 
The resulting integrals are substituted in Eqs. (IA1 U -( lA17b 



with the help of Eq. (lA24l - llA27l The final result is 
sented as 



P^(k) = exp {- [l + /(/ + 2)p 2 ] (k/k NL ) 2 } 

X (l + (F') + fp 2 f P h (k) + J] S'TEnmik) 



where 



and 



6;r 2 / 



dkP h (k) 



-1/2 
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(A69) 
(A70) 



9 „ 3^ 1 10 6 
+ <^">(^<25 + 2e 7 + ^i + yi?2 



+ <F"> (~ 68 + Qg) + {F'f [to + Qn + *Ri + 
1 



+ 2{F'){F")Q n + ^{F") 2 Q n , 



18 12 40 24 

En = — Qi + —Qi + 2Qi + — Rj + —R 2 

, /18 36 



(A71) 



+ (F")\-Q S +2Q 9 



+ a-'^\2Q 9 + 2Q u + ^R l + ^R 2 



2(F')(F")Q 12 , 



(A72) 



3 3 16 ,16 
En = - — Qi ~-Q 2 + -Qa- -Ri + (F } \Q 6 - -Ri 



\(F")Qs-^(F') 2 (Qs-Qio), 



(A73) 



57 51 1 16 30 

£22 =— fit + j^Qi + 3fi 3 - + yfli + y/?2 

, /12 18 24 

+ <^'> (yfis - Qe + 6g 7 + —Ry + —R 2 

+ (F")l~Q$ + Q 9 



+ {F'Y I - Q& + g 9 - -Q10 + Gn 



£23 = ~Q\ - 3Qi + \Qa - ^1 + (F')Q 6 , 



E24 ~l6 Qu 



(A74) 



(A75) 



(A76) 
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3 27 1 6 12 

£33 =jQi + —Qi + 2g 3 - + ijRi + —Ri 

+ {F')(-Q 6 +2Q-j), (A77) 
£ 34 = -|fii - jfi 2 + (A78) 

E^ = ^Q 1 + ^Q 2 + ^Q 3 -^Qa, (A79) 

and all the other E„ m which are not listed above are zero. 

Equations dA69HA79l ), together with Eqs. ( |A39l )-( lA48b . 
are complete set of equations to give the general one-loop 



power spectrum with effects of local Lagrangian bias and 
redshift-space distortions. Although the number of terms are 
large, they are all given by simple integrals of Q„(k) and R n (k) 
of Eqs. dA39l ) and dA40t . which are numerically easy to eval- 
uate. The power spectrum in real space is obtained by simply 
putting / = 0. When the bias is not present, (F') = (F") = 0, 
this result exactly agrees with the one which is derived in pa- 
per I. A cross power spectrum is obtained by substitutions of 
Eqs. dA181 l— dA22b . after expanding Eq. dA691 > in terms of (F'} 
and (F"). 
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